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Abstract. We compute the asymptotic expectation of the number of 
open nodal lines for random waves on smooth planar domains. We 
find that for both the long energy window [0, A], and the short one 
[A, A + 1], the expected number of open nodal lines is proportional to A, 
asymptotically as A — > oo. Our results are consistent with the predictions 
of Blum, Gnutzmann and Smilansky IBGS: in the physics literature. 



1. Introduction 
Let Q C M? be a smooth planar domain and consider the Dirichlet problem 

(1) - Aipj = \jipj; ifj\dQ = 0. 

Here, we assume that f^ipiipjdx = 5j. Consider the zero set of the j-th 
Dirichlet eigenfunction (with the boundary excluded) 

Nj := {x G O; <pj(x) = 0} - dQ. 

We call the set Nj the nodal set of the eigenfunction (ftj. This set is a curve 
which in general has self intersections. However, for generic domains [TJ] . 
a nodal set is a union of connected components consisting of closed loops 
homeomorphic to circles and open nodal lines homeomorphic to open inter- 
vals. In this paper, we compute the average number of open nodal lines for 
a random linear combination of Dirichlet eigenfunctions in various spectral 
intervals. This problem is generically equivalent to counting the zeros of the 
boundary traces of the eigenfunctions. The latter problem has obvious sim- 
ilarities to measuring the length of the interior nodal line, and our results 
show that the order of magnitude is the same. We recall that S. T. Yau 
conjectured that in all dimensions, the hypersurface volume should satisfy 
c\j < TL n ~ 1 {Nj) < CXj for some positive constants c, C depending only on 
(M,g) |Yll IY2j . The lower bound was proved in dimension two for smooth 
domains by Briming-Gomes [BGJ and both the upper and lower bounds 
were proved in all dimensions for analytic (M, g) by Donnelly- Fefferman 
[DF1 IDF2] . Recently, there have been important advances in understanding 
other aspects of the geometry of nodal sets for large eigenvalues, including 
asymptotic results for the expected number of nodal domains for random 
spherical harmonics |NSj . as well as for the distribution of nodal sets on 
general Riemannian manifolds without boundary [Z2J. 
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To state our results, we define the random linear combination of eigenfunc- 
tions (i.e. random waves) corresponding to the long range energy window 

(2) $a(*;A):= a m (x), 

A 3 G[0,A] 

and the random combination corresponding to the short range energy win- 
dow 

(3) Yl w( x )> 

where in both cases, a,j are (0, 1) Gaussian i.i.d defined on the sample space 

(4) N L (\):=#{\ k €lO,\}}; N S (X) := #{A„ £ [A, A + 1]}. 

-L S 

The nodal set of the random wave $ a ' is by definition the curve 

A/- a L < 5 (A) = {xen; <S> L a ' s {x- A) = 0} - 90. 

We are interested here in computing the asymptotics of the number of in- 
tersections of the nodal set with the boundary, d£l. Let 

T a L ' 5 (A) = card(Ar a i ' 5 (A)n^). 

Since for generic domains, a nodal set is a union of closed loops homeomor- 
phic to circles and open nodal lines homeomorphic to open intervals, the 
intersection number Ia' S (A), is therefore, almost surely, equal to twice the 
number of open nodal lines of the nodal set, M a ' (A); with probability 1, 
there are no multiple intersections with the boundary (see Lemma l3.ip . 

Given the random variables, X a ' (A), we define the corresponding expec- 
tations 

Z L ' S (X) = E1^ S (X). 

Our main result is the computation of the leading asymptotics of Z L,S (X) 
as A — » oo. For the following theorem, we need to impose a generic non- 
recurrence condition. We say that a point q G d£l is non-recurrent if the 
measure of loops at q G <90 for the associated billiard map j3 : B*dVl — > B*dVt 
is zero. In the following, we denote the arclength of dSl by 1(80,). 

Theorem 1.1. Let <fj',j = 1,2,... be the Dirichlet eigenf unctions of a 
smooth domain 0, and asssume that all boundary points are non-recurrent. 
Then, given aj;j = 1, . . . , N L ' S (X), i.i.d. (0, \)-Gaussian random variables, 

(<) Z L (X) = e -^Bx + o(X), 



and 



as X — > oo. 



(ii) ^(A) = ^A + o(A) 



Remark 1.2. The non-recurrence assumption in Theor em l 1 . 1 1 holds for generic 
boundaries and, in particular, for all convex analytic domains [Zlj . 
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The result in Theorem 11.11 is related to well-known results of Berard 
[Bej on the expected nodal lengths of random superpositions of eigenfunc- 
tions and Berry [BR], who computed the expected length of nodal lines for 
isotropic, monochromatic random waves in the plane (eigenfunctions of the 
free Laplacian). He also gave a somewhat heuristic argument for comput- 
ing the asymptotics of the variance. Of more direct relevance is the recent 
result by Zelditch [Z2j on the expected nodal distribution of random waves 
on compact manifolds without boundary. One can naturally view Theorem 
1.11 as the analogue for domains with boundary. Blum, Gnutzmann and 



og 

Smilansky [BGSjj have studied the distribution of the number of boundary 
intersections Vj of the nodal set of (pj (in addition to the number of nodal 
domains). Using Berry's random wave to model <pj for chaotic systems, they 
found that for large eigenvalues, Vj should concentrate around S- ■ Xj, where 
t is the length of the boundary dVt. Numerical results for eigenfunctions of 
both Sinai and stadium billiards support this prediction. Since for general 
domains, the average nodal asymptotics over spectral intervals [A, A + 1] 
should be the same as for individual eigenfunctions of ergodic billiards, our 
asymptotic result for Z s {\) in Theorem 11.11 is thus consistent with [BGSJ. 

Recently, in the case of piecewise-analytic domains, Toth and Zelditch 
[TZlJ have proved deterministic upper bounds for the number of nodal in- 
tersections with the boundary dfl (and more general interior curves) for 
individual eigenfunctions. In work in progress, when f2 is an ergodic bil- 
liard, these authors have also proved some asymptotic results for the nodal 
(and critical point) distributions of complexified restrictions of Dirichlet and 
Neumann eigenfunctions along strictly convex, real-analytic interior curves 
C C £1. In this case, at least for an ergodic sequence of eigenfunctions, the 
number of complex zeros of the holomorphic continuations of the eigenfunc- 
tion restrictions, (j)j\c, is ~ cA. The same result is likely true when C = d£l 
and this would of course be consistent with the random result in Theorem 

O 

In examples like the torus or sphere, where the spectrum of f2 is degen- 
erate with high multiplicity, rather than summing eigenfunctions belonging 
to different eigenspaces, it is natural to consider the ensemble of random 
eigenfunctions attached to fixed eigenspaces. A natural way to do so is to 
fix a basis B = {rji, . . . rjj^} of the eigenspace £\ and consider the random 
ensemble of functions on f2 defined by 

7] = aiT}l + . . . + (WT/Af, 

where dj are standard Gaussian i.i.d. Note that the probability density of n 
is independent of the choice of the basis B. 

Berard |Bej computed the expected length for the nodal line of a random 
eigenfunction on the sphere to be const ■ A. Rudnick and Wigman [RW] and 
Wigman |Wj have studied the variance of the length of the nodal line of ran- 
dom eigenfunctions with A — > oo, for the torus and the sphere respectively. 

Recently, Granville and Wigman |GW] have determined the asymptotics 
of the variance of number of zeros of random trigonometric polynomials of 
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degree ~ A and moreover, proved a central limit theorem for their distri- 
bution. While there are clearly similarities between the boundary traces of 
([2]) and ([3]) and random trigonometric polynomials of degree ~ A, it would 
likely be difficult to prove a central limit theorem for nodal distributions of 
random waves on arbitrary smooth domains. However, we do hope to study 
the variance and higher moments in future work. 

We thank Zeev Rudnick for many helpful discussions about random zeros, 
Steve Zelditch for helpful comments regarding pointwise Weyl laws at the 
boundary and the anonymous referee for helpful comments and suggestions. 
The second author would also like to thank the CRM analysis laboratory 
and its members for their support. 

2. A PRELIMINARY LEMMA ON NODAL INTERSECTIONS WITH 80, 

As above, we let n be a smooth bounded domain in M? and let ipj be the 
I? normalized eigenfunction of the Laplacian with eigenvalue Xj satisfying 
the Dirichlet boundary conditions. Without loss of generality, we assume 
that the eigenfunctions are real-valued and we let 

q : [0,£] ^ On, 6 ^ q(9) = ( qi (e),q 2 (e)) 

be the arclength parametrization of the boundary d£l with I := £{dVL). 
Let Vj : [0,1] -> 1 be the boundary trace of (fj, that is, 

(5) Vj (e) := d^M^))- 

Here, v = v{q) denotes the unit outward-pointing normal to the bound- 
ary. Recall that the main object of our interest are the real valued random 
variables ^(x;X) and 3>f(x;A) (see © and ©). 

We would like to compute the leading asymptotics of the expectation 
Z L ' S (\) = Kla ' (A) and we do this by counting the zeros of the boundary 
traces of and <I>f defined by 

(6) V a L (d;\)= Y ¥iW. 
and 

(7) V a s (9;X)= Y a i v M- 

AjS[A,A+l] 

The functions V a ' A) are useful, since, as we show in the next section, 
their zeros correspond to the intersections of the nodal sets M a ' (A) with 
the boundary <9f2, for a generic choice of a. In particular, we will show that 
the number Ia' S (X) of the intersections of the nodal line Na' (A) with the 

L S 

boundary equals the number of the zeros of V a ' {0;X), almost surely (see 
Lemma l2.ip . This observation implies that 

(8) Z L ' S (X):=E[l^ s (X)]=E[lt S (X)], 
where 

3-a' S {X) = #{# G [0,1] : V a L ' S (0;X)=O}. 
Our main interest here is in the asymptotic behaviour of Z L,S (X), as A — ► 
oo. In both cases we reduce the problem of counting the nodal intersections 
to counting the number of zeros of random functions, for which we use the 
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Kac-Rice formula. There are related formulas that have been used by several 
authors in different settings: In a series of papers starting from [SZJ such 
formulas are used to study the distribution of zeros of random holomorphic 
sections of vector bundles. A more classical reference for such formulas is 

2.1. Open nodal lines and boundary critical points. In the Dirichlet 

case treated here, the key to computing the A — > oo asymptotics of Z(X) 
is the observation that generically the number of boundary critical points 
of & a (9;\) equals twice the number of open nodal lines. The following 
elementary result makes this correspondence more precise. 

We begin with the following result for functions vanishing on d£l. It is 
an elementary result which is probably known (see e.g. |BGS] ). but since we 
could not find a direct reference, we include the proof here. 

Lemma 2.1. Let $ G C°°(f2) satisfy $>(q) = for all q G <9fL Assume that 
the critical points of $ on dfl are simple. Then the number of open nodal 
lines of <I> intersecting dQ equals \ the number of boundary critical points. 

Proof. Let A/$ be the nodal set of First, we assume that q is an intersec- 
tion of A/$ with d£l. To show that d v <&(q) = 0, we flatten out the boundary 
and introduce local normal coordinates {x\,X2) G (— e,e) 2 with x(q) = (0,0) 
and also, X2(p) = for all p G <9S1 with X2 > in 0,. 
We locally have either 

(9) (f(x 2 ),x 2 ) GA/i 

for some / G C°°(— e, e) with /(0) = 0, or alternatively, 

(10) (xi,g(xi)) G A/$, 

for a C°°(— e, e) function g with g(0) = g'(0) = 0, where g doesn't vanish 
identically on any interval containing (the latter case is if the nodal line is 
tangent to the boundary; only the former case is possible, see the proof of 
the converse statement below). In fact, the simple zeros assumption ensures 
that locally there is an equality rather than an inclusion (see the proof of 
the converse statement). 

In case ([9]), the Taylor expansion of &(xi,x 2 ) around (x\,0) gives 

$(f(x 2 ), x 2 ) = d X2 $(f(x 2 ),0)x 2 + 0{x\) = 0, 

since $(xi,0) = 0, so that 

d X2 $(f(x 2 ),0)=O(x 2 ). 

Since /(0) = 0, it follows that d X2 Q(0,0) = 0, which finishes the proof in 
that case. 

In case (jlOD . similar reasoning gives 

d X2 $(x 1 ,0) = O(g(x 1 )), 

provided that g{x\) does not vanish identically. Since g(0) = 0, it follows 
again that d X2 &(0, 0) = 0. 

To prove the converse, assume that q G d£l is a critical point of $ and 
show that q is a nodal intersection point with the boundary, d£l. Using the 



6 JOHN A. TOTH AND IGOR WIGMAN 

same normal coordinates as in the previous proof, we consider the equation 

(11) *(X!,X 2 ) = 0, 

where, by assumption d X2 &(0, 0) = 0. The fact that <£(xi,0) = implies 
that 

$(xi,x 2 ) = x 2 ip(x 1 ,x 2 ) 

with tjj (0,0) = and d Xl t/j(0,Q) ^ 0. It follows by the implicit function 
theorem that ip(xi,x 2 ) = is locally the graph x 2 = g{x\). 

□ 

3. A Kac-Rice formula for Z l ' s (X) 

3.1. Some notation. Given a point 9 G [0,^] and the spectral parameter 
A G M+, we introduce the vectors b^ S (9; A), b 2 ' S (9; A) G M Ari,S ( A ) where 

(12) bf(9; A) = (v k (9)) Xk€[0jX] ; b L 2 (9- A) = (d e v k (9)) Xhe[0>X] 
for the long spectral range, and 

(13) bf(9;X) = (v k (9)) Xke[x>x+1] ; b s 2 (9; A) = (d e v k (9)) Xke[x>x+1] 

for the short range. Here the dimension is given by (JH). Note that for A 

L S L S 

sufficiently large, for every 9 G [0, t\ , the vectors b 1 ' (9; A) and b 2 ' {9; A) are 
not collinear (see Corollary 16. ip . We will use this fact in the proof of the 
main result of this section (Proposition 13. 2| ) . 

In addition to the vectors b L,s {9; A) of boundary traces of eigenfunctions 

L S 

it is also useful to define the corresponding functions c-' (9; A) : [0, £] — > M 
for i,j G {1,2}, where 

c£(0; A) := (bf(e;\), bf(9;X)); cfj(9; A) := <&f(0;A), bj(9;X)). 

For example, c^f is just the squared length of b±' S and c^ 2 ,S = c 2 { ■ 

3.2. A Kac-Rice formula. In this section we give a Kac-Rice type formula 
for computing the expected value of the number of zeros of a random com- 
bination of Dirichlet eigenfunctions. Results similar to the one here can be 
found in [CL] (see pg. 285). Bleher, Shiffman and Zelditch [BSZlj . [BSZ2] 
have proved a Kac-Rice formula which applies in a more general situation. 
However, for the convenience of the reader, we give a direct elementary proof 
of Proposition 13.21 in the appendix. 

The following Lemma f3,ll is used in the proof of Proposition l3.21 It implies 
that with probability 1 all the zeros of the boundary trace are simple and 
together with Lemma ET] this implies the equality (JEj), i.e. that the expected 
number of nodal intersection with the boundary equals the expected number 
of zeros of the boundary trace. The proof of Lemma 13. II is also given in the 
appendix. 

Lemma 3.1. For A sufficiently large, the set 

C = {a G R AfL,s ( A ) : 39 G [0,4 Vf< s (9; A) = d 6 V a L > s (9; A) = 0} 
satisfies fi(C) = where dfi(a) := (2ir)~ NL ' S : ( A )/ 2 e~l' a ll 2 l 2 da. Moreover, 

coalimC > 1. 



COUNTING NODAL LINES 



7 



Proposition 3.2. The expected number of nodal intersections with the bound- 
ary d£l is given by 



4f(0;A) fc^(9;\)^ 2 



de, 



(14) Z L ' S (\) = 

o * 

L S i 

with (6; A) defined in section [3A[ 

In the following sections H] and we show that the leading term on the 
RHS of (fT4"j) is ~ C^ S A for some universal constants Cq S > 0. 

4. ASYMPTOTICS FOR Z L (\). 

In this section, we compute the large-A asymptotics of the pointwise sums 
Cii(0;A) = Ea,<>#)| 2 and 4(0; A) = £ Aj < A |^-(0)| 2 . The leading 
asymptotics for the special case of c±±(9; A) was computed by Ozawa |Oz] 
and the asymptotics for the pointwise sum C2 2 (#; A) is closely-related, but 
we could not find it in the literature. For the benefit of the reader, we 
give the computation of both of the diagonal sums here. Alternatively, the 
leading asymptotics for c^{9; A) and C2 2 (#;A) follow from wave analysis in 
section [5j However, since it is more elementary to use the heat calculus for 
manifolds with boundary (and also provides an independent verification of 
the asymptotics), this is the approach we take in this section. The required 
bound for the mixed sum cf 2 (#; A) is more subtle and requires the full two- 
term asymptotics arising from the wave analysis in section [5] (see (|49p ). 

The relevant heat parametrix construction goes back to work of R. Seeley 
[S] and here we briefly recall the main results that will be needed for our 
computations. The reader can find more detailed treatments in [HZ] Ap- 
pendix 12, [S] and |Uzj . Roughly speaking, the Seeley parametrix for the 
Dirichlet (or Neumann) resolvent for an elliptic boundary value problem is 
constructed as a sum of an interior parametrix and a Poisson- kernel- type 
correction which compensates for the boundary conditions. To describe it 
in more detail, it is useful to introduce normal coordinates (x, y) G R 2 in a 
neighbourhood of the boundary so that the Euclidean metric takes the form 
dx 2 + h(x, y)dy 2 and the domain is given in these local coordinates by x > 0, 
the boundary corresponding to x = 0. Consequently, since q : [0, £] — > dfl is 
the arclength parametrization, 

(15) d0 2 = h(0,y)dy 2 . 

We let (x, y) G U where U is a sufficiently thin tubular neighbourhood of 
the boundary, <9f2. 

In local coordinates, the N-th order Seeley parametrix for the Dirichlet 
resolvent Rn(i^) := (— — (here, \i = A 2 ) is of the form 
(16) 

N 

Rn,N(^) = E( 2vr )~ 2 / / e^-^e^-y'^c. 2 . :i (x,y,tv^)h- 1/2 (x,y)d^dr 1 

N 

+ V(2vr)- 2 / / e- ix ^<y-y^d^ j (x,y,(,r,, t j,)h- 1 / 2 (x,y)^dr 1 . 

~q JM 2 Jm. 2 
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In (p~6|) the factor h~ l l 2 {x,y) is included in the integrand to cancel the 
coefficient in the volume measure /i 1//2 (x, y)\dxdy\ upon integration. The 
first term on the RHS in (|16p is the interior parametrix and the second one 
is the boundary-compensating parametrix. Here, we are mainly interested in 
the leading coefficients, c_2,d-2- A direct computaton in local coordinates 
shows that 

(17) c_ 2 (x,y,Z,r J ,\) = (e + \v\l-Vr 1 

where, 

\V\1 : = h^ix^y)^ 2 . 
For the subsequent terms in the asymptotic series J2f=oC~2~j, 

sup iD^c^-jix^^^;-^)] < C a (\C\\ + \v\ + \^\ 1/2 )~ 2 ~ j ; J>1,H>0. 

Here, we abuse notation somewhat and denote the holomorphic continuation 
of C-2-j to the cone S = {(Cm) G C 2 x C;3f/i < + |^C| 2 - I^CI 2 )} 

with C(C) 7 ?) a l so by c -2-j- For the boundary-compensating terms one has 



e -*y/\e\l-v I — - — 

d-2[x,y,i,r],n) = - 2 , Sftv/ r?^ - /x > 0, x > 0. 

In this case, the subsequent terms in the sum J2f =0 d-2-j satisfy estimates 
of the form 

sup \x a DPd^ j {x,y^r ] -i l £)\ < C a ^\H\ + |r/| + ^V^-i-i-M+lfl 
x,y£U 

x exp(-<5 a)/3 y(|£| + |r?| + |m| 1/2 ) ); > 0. 

To get a parametrix for the Dirichlet heat kernel of £1, one writes the heat 
kernel as a contour integral 

(18) e-* A « = / e-*"(-A n - m) -1 ^- 

Here, T = {— L + se ±ie ,s > 0} with L > and 9 E (0, §), is a wedge 
enclosing the spectrum of A^. 

Substitution of the resolvent parametrix (|16p in the contour integral (|18p 
gives the following expression for the heat kernel parametrix in U x U : 
(19) 

N 

H N (t)(x,y;x',y')=J2(^r 2 / / S^^-if)^.^ y, £, rj, i) hT^x, y) d^drj 
i=o 7 J 

N 

+ ^(2tt)- 2 / / e-^e^-y^5 j (x,y,^r ] ,t)h~ 1 / 2 (x,y)d^d7 1 . 
j=o J J 

It is well-known that Hj^(t) is a good approximation to e~ tAn in the sense 

that for t > 0, 

(20) 

|^^(e- tA «(z,ii;)-^(t)(z,u;))| < Ca ^ n+a+ ^ +N+1 ^ 2 e- 8lz ~ wl2/t ; 5>0 
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where n = dimfi = 2. By the Cauchy integral formula, 

lo{x,y,£,v,t) = [ e- ttl c-2{x,y,t,ri,n)dn = [ & 2 dfi = e~ t{ - e+ ^h) . 

Similar ily, 

if if e - x \P^F^ 

= _ e -Ke+\v\i) x e ~ix$ 

So, subsititution of the formulas for 70 and <5o and differentiation under the 
integral sign in (fT9|) gives 

(21) Yl^^My))? = dxd x ,H(t)(x,y;x' ,y')\ x=x/=0 , y=y , + 0(t-l). 

3 

= (2vr)- 2 / l - 1 //2 (0)y) J J e -*S 2 ^^4 e -^l 2 , e -^5( e ^_ e -^)]| 2=K , =0 ^ dr7+O (^|) 

= (2vr)- 2 / l - 1 / 2 (0,y) | I d x [e~« {jX-e**) ]|, =0 ^+0(H ) 

In the last line, differentiation of e - *'''^ gives an integrand that is odd in £ 
and so it integrates to zero. Thus, 

^e- a ?|^^(0,y)| 2 = 2(27r)- 2 / l - 1 /2( ,y) / e ~« 2 £ 2 «X / e^l^dv+O^H) 

= a(y)t- 2 + 0(t-f), 
where a(y) = 7-. By the Karamata Tauberian theorem, 

(22) c&(0;A)= E l«iWI 2 ~^^i(e)M a = ^WA 4 , 
where, 

(23) c f lW = (8vr)- 1 . 

Similarily, the asymptotics for the second pointwise Weyl sum ^a 2 </^ \^eVj(6(y))\ 2 
follows from the formula 
(24) 

£ e~* A ^#)| 2 = (2^- 2 / l - 1/2 (0,y)|9 e y| 2 4^^^ii'(t)(x,y;x , ,y / )]U-=x=o,j / = y ' 

= 2(2vr)- 2 / l - 1 /2 ( o, y) |^ y | 2 | | e-^^^f^lOtr 5 / 2 ) ~^ 0+ &(j,)r 8 , 

where, b(y) = ^^ 1/2 (0, y)/i(0, y) 3 / 2 1 <9 e2 /| 2 . 

There are Jacobian factors |<9gy| > occuring in (|24p because of the 
change of parametrization of the boundary given by y 1— > 6(y) and we also use 
that dgx = 0. Differentiation under the integral sign of the heat parametrix 
is justified since all exponential sums of derivatives of boundary traces of 
eigenfunctions are absolutely convergent in view of (|20p and derivatives of 
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the remainder term are controlled. Just as in (|22p . the Tauberian theorem 
gives 

(25) 4(0; A) = J2 \d0Vj(O)\ 2 ~A-,oo 4(#)A 6 , 

A,<A 

where, 

(26) c&(0) = ^[h(0,y(9))] \d e y\ 2 = (487r)-\ 

In ([26]) we have used that 9{y) = f£[h(Q, s)] 1 / 2 ^ (see (|I5]> ) and so, |%y| = 

[Mo,^))]- 1 / 2 . 

Finally, we turn to the mixed sum Ci 2 {9; A). By the same heat analysis as 
for c 22 (9; A), one shows that 

\devj(9) + vj(9)\ 2 -a^oo (48tt)- 1 A 6 

and so, by writing 2cf 2 (0; A) = £ Aj < A \iPo + 1)«;0?)| 2 - Ea 3 <a l^i W ~ 
T,\ J <\\vj{6)\ 2 it follows that cf 2 (0;A) = o(A 6 ) uniformly for 9 6 [0,4 Un- 
fortunately, this bound is not sufficient to prove our theorem. Indeed, we 
will need the stronger bound 

(27) cf 2 (#;A) = o(A 5 ). 

This estimate follows from analysis of the Dirichlet wave kernel and we defer 
the proof to section [5] (see (T39j) ), 

5. ASYMPTOTICS FOR Z S (\) 

5.1. Computation of cfi(9; A). The heat analysis in the previous section is 
not sufficient to deal with the short-range case, even for the diagonal sums. 
An alternative approach involves the Dirichlet wave operator. However, un- 
like the boundaryless case, one does not have an explicit wave parametrix, 
even for small time. This is the main complication in the case of manifolds 
with boundary. Indeed, the small-time behaviour of the Dirichlet wave op- 
erator kernel En(x,y;t); x,y £ is quite complicated and in particular, 
has no conormal expansion at t = 0. Nevertheless, it turns out that the 
restriction to the boundary diagonal does have such an expansion (see ([30]) ). 
Using the fact that the singularity of E^(q,q;t); q G d£l at t = is iso- 
lated, Zelditch |Zlj computed the asymptotics of cf^(9; A). The idea is to 
use boundary wave front calculus and the conormal expansion of Ivrii [I] 
and Melrose [MlJ together with the Hadamard variational formulas for the 
Dirichlet (or Neumann) eigenvalues to compute the coefficients in the conor- 
mal expansion of the boundary trace of the interior wave group restricted 
to the boundary diagonal. The asymptotic expansion of c 22 (9;X) and the 
bound for cf 2 (9;\) in the next subsection are apparently new. Since the 
computations for c 22 (#;A) and cf 2 (0;A) are related to the one for cf x (0;A) 
(and use it), we give the argument for each of spectral sums. 
Let 

E n (t) = costy^o 
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be the Dirichlet wave operator where denotes the Dirichlet Laplacian on 
O. To fix notation, we denote the boundary restriction operator by 7,9^1 s ° 
that for u G P'(fi) with WF(u) n iV*<90 = 0, y m u is the uniquely-defined 
restriction (ie. pull-back of u) with jgnu G 2?'(<90) and 7gn?i = ^|an when 
u G C°°(0). 

The diagonal kernel t n En(t, x, x) ^ C°°((0, e] x x $7) for any n > and 
e > when 2; G 0. The singularity at t = is not even isolated due to 
reflecting loops based at x G with period t = 2d(x,dQ). However, when 
x = q G 90, the singularity at t = is isolated. The reason for this can be 
described as follows: let 

$ t : T*0 -> T*0 

be the time-t broken generalized bicharacteristic flow in O MSj . For (q, 77) G 
£?*<90 we let £(q,T]) G S+(0) be the unique inward pointing unit vector 
that projects tangentially onto 77. We denote the tangential projection of 
v G S+(0) here by v T G B*(dQ) and 79^ : <p l— > ^0|an the boundary trace 
operation. 

Propagation of singularities in O |MSj[ IM1| gives 

(28) WF(Ea(x, y, t)) = {(t, r, x, £, y, rf)\ t = \£\ g , $ t (x, = {v, v)}- 

For fixed t G R we define the distribution Eqj on Q x 11 given by 

En,t(x,y) = E n (x,y,t). 

Let N*Aqq = {(q,£,q,—£);q G <90} be the conormal space of the diagonal 
A-on C (90 x <90. In view of (f28|) . when i G (0, e) with e > sufficiently 
small, 

WF(E n , t )riN*(A dn ) = (/). 

This just says that for \t\ sufficiently small, there are no non-trivial, periodic 
broken bicharacteristics passing through q G <90 with period \t\. Then by 
wave- front calculus, E^ := jAdttEn is well-defined for e > small and for 
t€(Q,e), 
(29) 

WF(E\q, q, t)) = {(t, r, q, r], q, rf) G T*(-e, e) x IT 90 x £*90; r = \ V \ g , 

[MQ,C(q,v))} T = (q',V% 
In (|29|) . £(q,r)) G S* n St is the inward-pointwing unit co- vector at q G 90 
with [£(g, r/)] T = 77 where t> t is the tangential projection at q G 90. 

It follows that the singularity at t = is classical conormal [U IMlj , so 
that for e > small as t — > + , one has t n E b (t, q, q) G C°° where n = dimO 
and so, 

00 

(30) t n E\t,q,q)~J2ak(q)t k 

k=0 

with dj G C°°(90). Under the non-recurrence assumption on the billiard 
flow, the expansion in ([30]) holds, modulo C°°((0,e -1 ) x 90 x 90), for all 
t G (0, e _1 ) where e > is fixed arbitrarily small. 

Since the method of proof in (TJ IMlj is non-constructive, one needs an 
additional argument to compute the actual coefficient functions ctj. One 
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way to do this is to integrate both sides in (|30|) against a test function 
tp £ C°°(dO). This is allowed because t n E^> £ C°° and the result is that 

(31) t"£VW l^0i(g)| 2 V(g)^(g)~^O+E* fc / a k (q)^(q)da(q). 

j J «* k=Q Jan 

Consider a normal variation of the domain Q with variation vector field, ipv. 
By taking the variation 5^ of both sides in (j31f) and using the Hadamard 
variation formula for the eigenvalues: 

S^X 2 = 2Xj6^Xj = / ip{q)\d v ^)j{q)\ 2 da{q), 
J an 

it follows that 



z <9f I zt / zt 

J 13 



(32) =Y,e^ t f ^Id^iqtfdaiq). 

j Jan 

By a well-known asymptotic expansion for the wave trace pj EH] 

(33) ^ ~^o+ (27r)- 2 z;oZ(^)(t + iO)' 2 + C'vol{dn)(t + iO)' 1 + ... 
j 

where the dots denote lower-order terms. Substitution of (f33l) in (|32|) gives 
(34) 

™ f|E^( e% *)j ~*-o+ 2(2+l)(2vr)- 2 ^z;o/(^)(t+z-0)- 4 +C / ^t;o/(^)(t+i0)- 5 + 

Equating coefficiens in (f3"T]) - (f3"3"|) implies that 
(35) 

Ve iA ^ / V(g)l^0j(g)| 2 da(g) ~ t ^ + 6(2ir)-%vol(Q;)(t^0)^+CfS^vol(dQ)(t-H0y 

where, a direct computation gives 

5^vol{£l) = / ^(q)da(q). 

Under the assumption that there are measure zero loops at the boundary it 
follows from (|31|) by a standard Tauberian argument (see for example |Zlj 
section 5.2) that 

(36) c s u (9,X)= Y, b#)| 2 = 4W-A 3 + o(A 2 ), 

and since the test function ip £ C°°(5$7) is arbitrary, it follows from ([35]) 
that 

(37) 4(0) = (27T)- 1 . 
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5.2. Computation of cf 2 (#;A). The computation of cf 2 uses the asymp- 
totics for cf x above together with some additional integration by parts and 



applications of the local Weyl law |HZj for boundary traces of eigenfunc- 
tions. Let dg q : C OD (d^l) — > C°°(d^l) be the tangential derivative given by 
de q ■= dq{d e ). Since WF(d 9q ) C {(<?,£;<?,-£); (?,£) G T*dn}, by wave front 
calculus, 
(38) 

WF([d eq End e/q ] b (q,q',t)) = WF([d dq El l de' q }(q,q , ,t)) C WF(E b n (q,q' ,t)). 

Then by the same argument as in the previous section, it follows that 
dg q EQd$ q has a unique restriction to the boundary diagonal which we denote 

by [d$ q Efid$i](q,q,t) := lA an [de q EQd$i]. Since wave fronts restrict, under 
the non-recurrence assumption, for t £ (0, e -1 ) with e > arbitrarily small, 
modulo C°°((0, e _1 ) x d£l x d£l), one has the following conormal expansion 

00 

(39) [d eq E b n d e ,](q,q,t) ~ t _ 0+ E b k (q)r n - 2+k 

k=0 

where b k £ C°°(<9r2) and the corresponding two-term asymptotic formula 
£ l^#)| 2 = 4(0)A 5 + O (A 4 ). 

A 3 <A 

The remainder of this subsection is devoted to computing the leading coef- 
ficient cf 2 (0)- 

From now on, we put Hj = = 1,2,3, ... and use semiclassical pseu- 
dodifferential calculus on dQ. Let = Op^(a) with a G S^(T*dQ) and 
:= q^ l oOph J {a)°q be the corresponding semiclassical pseudodifferential 
operator on the parametrizing circle R/^Z. It follows that in the Dirichlet 
case considered here [HZJ, 



Aj<A 

where, iVo(A) := (2n)~ 2 (volB*Q,) A 2 . Since we are assuming here that the 
measure of periodic broken bicharacteristics is zero, using the conormal ex- 
pansion in (f39l) one gets that for some constant C a (yet to be determined), 

(41) £ ^<4,«W = C« + olA- 1 ). 

Aj<A 

But then from the weak law in (|4U|) it follows that C a = J b *qq a(y, r])da(y, rj) 



where da(y, rj) := 2(ir volfi) 1 y 1 — |r/|^ ch/dr/. Taking sums with Xj < A + 1 

and Aj < A in (|41|) and subtracting them gives 

(42) 



3 ' <<j 
Aje[A,A+l] 



hj{A% Vj ,Vj) = / a(y,r ] )da(y,r ] )[N (X + l)-N (X)] + o A (X). 



B-an 
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Now rescale and write Xj = hj 1 and let ip G C°°(K/^Z). In analogy with 
the previous section, we want to compute 

(43) Y f m\dev 3 {e)\ 2 d6= Y, K?ty{KjD 9 )\,Vj) 

+ Y hj 2 {(hjDQil>)vj,hjDevj) 
= Y ^(#^.^ + 5 E <^M 2 > 

The second last line in ([43 j) follows by integration by parts and the last line 
by the result in the previous section for 0^(0; A) (here, we also use that 
hj 1 G [A,A + 1]). Since 

WF' n .{F{h- 1 )) C B* 1+5 (dn) x B* 1+s (dQ), 

from the boundary jump equation (I57p it follows that for any 5 > 0, 

^.(^)c(^r 1 B* 1+5 (dn). 

Let X G C co ((^)- 1 J B 1 * +2(5 ^) with X (0,0 = 1 for (9,0 G (dq^B^dil. 
Then, from the last line in (|43|) . 

aj^iA.A+i] 

£ ^7 2 (ch (x)^(^^) 2 OpRj (xK,«i> + o(A 3 ) 

= A 4 (l + ©(A" 1 )) Y tfiOPKiixMfrjDefOpHjMvj^j) + 0(A 3 ) 

RT 1 6[A,A+l] 

Clearly, il>Op h] (x)(h j Dg) 2 Op hj (x) G Opn^S^) has principal symbol 

^[^opR,(x)i^l 2 PPB,(x)](f,0 = v>wi£l 2 , 

for (9, G (dq t )- 1 B*d£l. But then, since /i" 1 G [A, A + l], it follows from (@2J) 
that one has the two-term asymptotic formula cf 2 (^; A) = C2 2 (#)A 5 + o(A 4 ), 
where, 
(44) 

cg,(0) = lim A -5 — ^— f / |e| 2 VlHli 1 * J A 4 [iV (A + 1) - iVo(A)], 

A^oo 7TV01SZ W|^|<1 / 

and since iV (A + 1) - N (X) = 2(2tt)- 2 (vol B*VL) A, it follows that 
(45) cf 2 (0) = 2 x - x — \- x (27r)" 2 volS*fi = (Svr)" 1 . 



COUNTING NODAL LINES 15 

5.3. Bound for cf 2 (9;X). By the same argument as for c 22 (9; X) one gets 
the two-term asymptotic formula, 

(46) ]T \d eVj (0) + v,{9)\ 2 = (Svr)-^ 5 + o(A 4 ), 

Aje[A,A+l] 

uniformly for 9 G [0, £]. Then, 
(47) 

2cf 2 (9;X)= Idevm+v^l 2 - £ \d eVj (9)\ 2 - £ \ Vj (9)\ 2 

AjG[A,A+l] A 3 e[A,A+l] A 3 e[A,A+l] 

where, the final bound in (|47p follows from the two-term asymptotic formulas 
in ([37]), gS]) and gTj). The upshot is that 

(48) cf 2 (#;A) = (A 4 ), 
and so, by the triangle inequality, 

(49) |cf 2 (0; A)| < C n X \cf 2 (8; A)| = o(A 5 ). 

6. Concluding the proof of theorem 11.11 
Corollary 6.1. There exists a constant Xq > such that for A > Ao the 

L S L S 

vectors b 1 ' (9; A) and b 2 ' {9; A) are linearly independent for all 9 G [0, £]. 

L S 

Proof. By Cauchy-Schwarz, the vectors b 1 ' (9;X) = (vi(9), v \(9)) and 

L S 

b 1 ' (9;X) = (dgvi(9), ....,dgv\(9)) are collinear for some 9 G [0,£] if and only 
if 

c^ S (9;X)-c^ S (9;X)-\c^ S (9;X)\ 2 = 0. 
But by our pointwise Weyl sum computations in (|23|) . ([26]) and (]2"7j) . 

cn(0; A) • cg(9; X) - \cf 2 (9; X)\ 2 C n (9)X 10 ; C n (9) > 0, 

and similarily, from ([36]), ([M]) and ([15]) . 

c s n (9; A) • cf 2 (0; X) - \cf 2 (9; X)\ 2 ~ X -+oc C Q (9)X S ; C n (9) > 0. 

So, for A > large, neither expression vanishes. □ 

Remark 6.2. We note that by a similar analysis to the one given in sec- 
tions H] and [5] and Corollary 16.11 above, one easily shows that the vectors 
(dgvj(9))x je [x,\+i] and (d 2 v j (8)) X:je [x,x+l} are linearly independent for each 
9 G [0,1]. The analogous result for the long-range case is also valid. 

Proof of theorem Combining the asymptotic formulas for c^- (9; X) , i,j = 
1, 2 in ([23D, ([26]) and ([27]) and ([37]), gS]) and gSJ) gives 



qf(M) 



o(A) 



uniformly for 9 G [0, £|. On the other hand, 



where, 



c 2 2 (M) c 2 2 (#) X 2 

cii (5; A) c x { (9) 



40) _ i and 4(0 _ i 
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So, substitution into the Kac-Rice formula in Proposition f3T 



rZ " 5 (A) = ~| / dB\ A + o(A), 




implies that Z b (X) = ^A + o(A) and Z L {\) = ^A + o(A). 

□ 



Appendix A. Proof of the Kac-Rice formula 

Proof of Lemma \3.1\ The proof in the long and short range cases is the 
same, so without loss of generality we treat the long range case here. 
In the following, we fix A > and put N = N L (X). Consider the map 

* A : R N x [0, i] -> R 2 
(a,9)^{V a L (9;\),d e V a L (8;\)). 

Clearly, 

C = vr 1 (^ 1 (0,0)), 

where tt\ is the projection onto R N . We claim that ty\ is a submersion. 
Given this, one has 

dimtf^O.O) = N - 1, 

and so, dimC < N — 1, which proves the Lemma. 

To see that is a submersion, we note that its (N + 1) x 2 differential 
matrix is given by 

where frf^A) and 62(^5 A) are the vectors introduced in section ETT1 The 
matrix is of full rank (i.e. rank 2) for A sufficiently large, by Corollary 
16.11 and so ^ \ is a submersion. 

□ 

Definition 1. Let X[— 1,1] be the characteristic function of the interval [—1, 1] 
and fix e > 0. We introduce the random variables 

(50) % S (A):=^/x Mll] (^^)|^(^;A)|^. 



L S 

When V a ' (•; A) has no double zeros, 

hm l^ s (\)=lt S (\). 

To compute ETq' (A), we need to interchange the e — > + limit and 
integration over a G K^' This requires showing that Ia j , f(X) is bounded 
uniformly in e > by a /x-integrable function in the o- variables. 
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~ L S 

Lemma A.l. Fix A > sufficiently large and let I a ,e (A) be defined as in 

H5DD- 

(i) In the case where dtl is C°°, 

«nn f L > S (\\ a T 1 fw NL ' S W- r -||a|| 2 /2 rf « A 

(ii) In the case where dVt is C u , one has the explicit bound 

sup X a L f (A) < C^ S X, 

where, C^' S > are constants depending only on $7. 

Proof, (i) The argument for both the long and short ranges are the same, so 
without loss of generality, we assume here that Xj E [A, A + 1]. To prove the 
first part of the Lemma, only the behaviour in a E M nS is important (and not 
in A). So, we henceforth fix A > sufficiently large and suppress dependence 
of all sets, functions, etc. on A. In the smooth case, the asymptotics in A 
of the number of boundary critical points is, to our knowledge, open and 
probably rather difficult (see however part (ii) and |TZlj for a sharp result 
in A in the real-analytic case). 

First, we note that since Vf(9;X) = Ea e[a,A+i] a j v j(@) * s linear in a £ 

~R nS , the number of boundary critical points, C(V a ;X), is invariant under 
scaling by ||a|| ^ 0; that is, 

(51) C(V a ;X)=C(V u ;X) 
where, u := tA-. Thus, we are reduced to proving that 

(52) C{V u \\)eL\S NS -\dtJ), 

where du is the standard, unit constant curvature hypersurface measure on 
§ NS -\ 

The second point is that by the same argument as in the previous Lemma 
13.11 (see also Corollary 16.11 and Remark [6. 2p . with V a (resp. dgV a ) replaced by 
dgVu (resp. dgV w ) one shows that given the map de^\ : S^ -1 x [0,£] — > R 2 
defined by 

d e ^x : (co,9) ^ (deV u (9;\),ti$V u (0',X)), 

the set dg^ 1 (0,0) C _1 x R is a finite union of compact, C°° hypersur- 
faces (here, we use that V u (0+£; A) = V LU (9; A) for all 9 E M). But then, since 
_1 x R = R w — 0, by the generalized Jordan-Brouwer separation theo- 
rem applied to each of the compact, connected hypersurfaces [2], it follows 
that (S N - 1 x R) - dg^ 1 (0,0) has finitely -many connected components. 
Now consider 

C := {u E S* 5 " 1 ; 39 E [0,£},d e V UJ (9; A) = d^V u (9; A) = 0}. 

Written another way, C = tt^^^O, 0)) where vr : S^" 1 xl-t S 7 ^ 1 is 
the smooth canonical projection map tt : (u>, 9) t— > u. Since ir maps connected 
sets to connected sets, it follows that 

fl°(S JvS - 1 -C / )<c» 
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and clearly, uj(C) = 0. Now, make the decomposition §> N ~ x —C = B%Li- ■ - U 
Bm] M < oo where the -B/s are the (open) connected components (which, 
implicitly depend on A.) Without any loss in generality, we choose a point 
uq £ B\. Then, all boundary critical points of V LJQ (6; A) are simple and we 
denote them by 80,..., 9 p where P is finite (again, implicitly dependent on 
A). Since the following argument is the same for each of the critical points, 
we consider the first one: 9 = 9q. For (0 O , ojq) £ [0,£] x S N -\ we have by 
definition 

(53) d e V wo (e ;\) = 0, 
and since ujq £ B\ C S nS ~ 1 — C , it follows that 

(54) c^ o (0 o ;A)^O. 

Then, by the implicit function theorem, it follows that for u £ B±, there 
is a unique C°° family of solutions 

(w,0oM) G B 1 x [0,£] 

to ([53]) satisfying 6>o(u;o) = #o- Repeating the same implicit function theorem 
argument for the other zeros gives the existence of smooth families of so- 
lutions 9i(u), 0p{oj) to (|53p for u £ Bi with respective initial conditions 
9±,...,9p. One can apply the same analysis to each of the other connected 
components B2,-.-,Bm- Let 

Bk, m = {lo £ B k : C(V U ; A) = m}, 

so that 

00 

m=0 

The argument above implies that each of the B^ m is open, and thus also 
closed (being the complement of |J Bfc m ') in B^. This implies that for 

m'^m 

each m, either i3fc,m = -Bfc or Bk jTn = 0. We conclude that C(V U ;X) is 
constant (and finite) on each of the (finitely-many) connected components 
B k of the N - 1 sphere. It follows that C(V U ;\) £ L 1 ^ 8 ' 1 ; du) and so 
also, 

So, given (i55|) . it suffices to prove that Xf e (A) = 0(C(F ;A)) uniformly 
in e > 0. The domain of integration in the definition (I50p of Xf (A) is a finite 

n 

union [0,£] = U Jj of disjoint intervals Ij = [(cj,dj)] each containing at 
3=1 

most one critical point. It is clear that the contribution of each of the Ij 
to the integral (|50p is at most 2. Moreover, for each j one of the following 
holds: 

(1) Cj = 0. 

(2) dj = £. 

(3) (0; A) = for some 9 £ Ij. 

(4) Either (V a s ( Cj ;\) = -e and Vf(dj;X) = e) or (Kf( Cj ,A) = e and 
y a 5 (^;A) = -e). 
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In the latter case, if j < n, then V£'(6;X) = for some 6 G [dj,Cj + i]. 
Therefore, n < 2 + C(V a ; A) + (C(V a ; A) + 1) = 2C(V a ; A) + 3. Thus the inte- 
gral in (|50|) is bounded by 

Zj e (A) <2-n<2-(2C(V a ;X) + 3) = 0(C(V a ;X)). 

(ii) Although we will not need the much stronger analytic bound in Lemma 
IA.1I (ii) in the current paper, we give the proof here since in light of recent 
work of Toth and Zelditch [TZ1| . we think it is of independent interest. 
From |TZlj Theorems 1-3, one has the following bound for the number of 
boundary critical points of individual eigenfunctions, Vj\ 

(56) C( Vj ) < C n Xj. 
Let F(A) : C7°°(dft) -» C°°(dSl) be defined by 

F(X)f(q) = [ d Uq G (q,q';X)f(q')da(q') 

where Go(q, q'; A) = IHa^^Ajg — q'\) is the free outgoing Greens function for 
the Helmholtz equation in M 2 . Roughly speaking, the bound in (I56D follows 
by holomorphically continuing both sides of the jumps-equation 

(57) d^jlan = -2F{Xj){d v ^\ m ), 

and using a Jensen-type argument to bound the number of complex (and 
hence real) zeros in a complex tube d£lc containing dtl by the exponential 
growth exponent of the holomorphically continued F(A)-kernel. Let Vj de- 
note the holomorphic continuation of Vj to complex parameter strip A{e) : 
[0,^] x [— e, e] where the holomorphic continuation of q is q c : A(e) — ► dQc. 
Then, one has the exponential growth estimate [TZlJ 

(58) sup |^ C (C)I < exp[Cf, |9Cl A^]. 

C6A(e) 

Moreover, 

V a C (C) = £ a jV f(C) 

Aje[A,A+l] 

and so by (I58p and Cauchy-Schwartz, 

(59) |K C (C)| < ^NS(X) \\a\\ exp[C n |9?C| A]. 
Since N S (X) ~ A, Theorem 3 in [TZlj applied to ^ gives 

(60) C(V a ; A) < C n max | log |V^(C)| \ < C' n (X + \ log ||a|| |) 

in view of (|59p . But by the scaling invariance (|5ip . it suffices to assume 
that ||o|| = 1, so the last bound in ([60]) is Oq(X). Now, just as in (i) one 
uses that sup e 2"^ e (A) = 0(C(V a ; A)). Again, the same argument works in 
the long-range case. 

□ 



Proof of Proposition \3.2[ We use (|8|) so that we are to compute the expected 
number EJ L - S (A) of the zeros of Vf' S (-, A), defined by © and © on [0,4 
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rL,S 



By Lemma [XT] we can assume that V a ' (9; A) has no double zeros so that 

V a L ' S (9;\y 
" — 

o 



\d,Vt s {6;\)\M, 



and so, 



E^ 5 (A) =Elimf / x\-i 



1 



Ye J X ^V 



o 



|^/" s (e;A)|de. 



Since by Lemma [A. 1} for each fixed A G M + sufficiently large, the function 

2 da 



supjff (A) G L 1 (iR^W; e-W 2 /2 



( 27r )^> s (A)/2 J ' 

by dominated convergence, we interchange the order of the integration and 
limit and get 

t 



ET a L ' 5 (A) = lim^-E J X [-i,i][ 



1 



fV a L ' S (0;X) 



d e V a L > s (9;X)\d9 



lim / E 





<J0, 



by Fubini. 

We rewrite the last equality as 



(61) 



EJ^(A) = lim / K^(0;X)d9, 



where 

(62) K^' s (9; A) = E 



2c 



Assuming A is fixed, we denote N = N L > S (X). To compute Kt ' (9; A) for 
a given 6 G [0,£], we note that 



b L l > S (9;X),a)=V a L > S (9;X) 



and 



(b^ S (9;X), a) = d e V a L >»(8;\) 



L,Sfn. 



By Corollary 16.11 for A large, b l ' (9; A) and b 2 ' (9;X) are nowhere collinear 

L S L S 

and so the vectors {6 1 ' (9; A), b 2 ' {9; A)} can be extended to a basis 

{b^ S (8; A), b^' S (0; A), b^ S (9; A), . . . , b L /(9; A)} 
of R N with the property that {b^ S (9;X), b^ S {0;X)} is an orthonor- 



mal basis of span{of' S (6>; A), b 2 '*(9; A)} ± . Let B L > S (9;X) G M N (R) be the 

L S 

matrix with row vectors b k ' (9; A). Then 



B L ' S (9;X)B L > s (8;Xy 



,,t_ (C L > S (9;X) 







In- 
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where 

rL,s ( a. \\ _ | c[f{9] A) c£f (0; An 
\c 2 { (9;X) c 22 (#;A)/ 
with Cij being defined in section 13.11 In particular, 



detB L ' S (9;X) = y det C L ' S (9; A). 

Writing the Gaussian probability density explicitly in (j62[) yields the for- 
mula 
(63) 

K e' (^; A) = — y X[_i,i] ^ JKV (6»; A), a)| exp (--||a|| ) ^ N/2 

We change the variables v = aB. In the new coordinates, we have 
||a|| 2 = a • a* = vB L ' S (9; Xy 1 (B L ' S (9; A)^ X )V 



( CL ' S ^ A) /jv °_ 2 ) = w\C L,s (9; \y l w\ + \\w 2 \ 



2 



where w\ = {v\,v 2 ) and w 2 = (t>3, • • • , f at), so that 
Kt S (9;X) ' 



(2vr) 7V /VdetC7 i ' 5 (e;A) 
J (y)!^! exp (- ^(wiC^^jA)- 1 ^ + ||i02|| 2 ))dt; 1 ---dvjv 



1 



(2vr) V / det C L > 5 (#;A) 

IK" 

/" / I,, ,|2\ du; 2 

/ exp(- -\\w 2 \\ )—n-2 ■ 



J ^X[-i,i]\^\v 2 \exp(-^w 1 C L ' s (9-,X)- 1 w t 1 )dwiX 



(2tt)- 



2 



Note that the last integrand is just the standard multivariate Gaussian prob- 
ability measure, so that the corresponding integral is just 1. 
Therefore, we get the formula 

(64) 

K^ S (9;X) = 1 =x 

27r v / detC L > 5 (0;A) 

00 e 

X YeJ /N ex P( _ \{v u v 2 )C L < s {e;X)- 1 (v l ,v 2 ) t )dv 1 dv 2 . 

— 00 — e 

We wish to apply the dominated convergence theorem again to exchange the 
limit and the integral in (|6ip . For this we estimate K^' (9;X) from above 
in the following way. Let E = {e^} be any orthonormal basis of Mr, with 
ei = p^jy- Then for a = Ea' the integral in (|64p is, using the invariance of 
the Gaussian measure, 

1 f f {b[' S {9, A), a^ei) ^ , , h L,S (a u P i\ 1 / 1 „ 
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We have 

(V (9;X),a 1 e 1 )=a 1 \\b 1 ' {6;\)\\ = a^^fc^, 
so that we integrate for \a'i\ < / T e a , =) an d we have by the Cauchy- 

Schwartz inequality 



(b^ b (9;X),Ea')\ < \\b%> b (0; \)\\ ■ \\Ea'\\ Vc 22 (fl; A) • ||o'| 



so that 



L.Si 



V c 22 (M) /■ 
— 2I — / dai 



1 



/„„-±||a'||2 



. i ii_/i|2 dai, . . . da'j 



e 4i 



N 



{2tt) n / 2 



\Rl 5 (0;A) : 



for some constant C\ > 0, since xe * x is bounded. Thus we have 
K^ ,s (6; A) < C\ 



(65) 



<C 2 A, 



for some constant C 2 > 0, where the last estimate in (|65p follows from the 
asymptotics for the c^'s in section 4 (in the long range case) and section 
5 (in the short range case). Thus by (|6ip and the dominated convergence 
theorem, 



E1 L ' S (X) = J K L ' S {6;\)d9, 



K L ' S (6;\) := lim K^ s {9;\). 

e— >0+ 



(66) 
where 



The fundamental theorem of the calculus and (JMJ) then imply that 
(67) 

CO 



-co 

oo 



2vr^det C L < S {6;\) 



1 c{f(0;A) 



1^2 1 exp(— - 



2detC L > s (0;A) 



-oo 

OO 



— OO 

l^(e;A)4 5 (^;A)-4 5 (0;A)2 



7T 



c£ S (0;A) 
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since / exp(— ^z 2 )\z\dz = 2. □ 
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